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An  experimental  and  theoretical  study  on  the  cathode  flow  pulsation  in  a  polymer  electrolyte  membrane 
(PEM)  fuel  cell  is  performed.  A  10-cell  PEM  fuel  cell  stack  with  open-air  cathode  channels  is  employed  to 
investigate  the  effects  of  the  cathode  flow  pulsation  on  the  overall  performance.  The  polarization  and 
corresponding  power  curves  obtained  show  that  both  the  limiting  current  density  and  the  maximum 
power  density  are  substantially  enhanced  when  the  pulsating  component  is  added  to  the  cathode 
mainstream  flow.  The  flow  pulsation  at  Re  =  77  provides  the  maximum  increment  of  40%  and  35.5%  in  the 
limiting  current  density  and  in  the  maximum  power  density,  respectively.  The  enhancement  of  the  overall 
performance  is  more  pronounced  at  low  Reynolds  numbers.  Also,  the  theoretical  mass  transport  analysis 
in  the  pulsating  cathode  flow  channel  is  carried  out  to  verify  the  present  experimental  results.  The 
momentum  and  species  conservation  equations  are  analytically  solved,  and  the  effective  time-averaged 
dispersion  coefficient  is  defined  to  account  for  the  enhanced  mass  transport  by  the  flow  pulsation. 
Comprehensive  analytical  solutions  show  that  the  effect  of  the  relevant  parameters  is  in  well  accordance 
with  the  experimental  results. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

As  a  new  and  environment-friendly  energy  source,  fuel  cells 
have  been  a  subject  of  intense  research  interest  due  to  their  high 
energy  efficiency,  low  emission,  few  moving  parts,  low  noise  and 
cleanness  of  energy  production,  to  name  a  few.  The  fundamental 
aspects  of  fuel  cells  have  been  well  documented  in  [1—3].  Among 
various  types  of  fuel  cells,  a  polymer  electrolyte  membrane  (PEM) 
fuel  cell  is  expected  as  one  of  the  most  promising  candidates  for 
future  power  source  on  account  of  its  high  power  density,  simple 
structure,  quick  start-up  and  easy  operation  [4]. 

A  practical  PEM  fuel  cell  system  consists  of  a  fuel  cell  stack  and 
a  balance-of-plant  (BOP)  system  for  the  stack  operation  such  as 
a  fuel  supply  subsystem,  an  air  supply  subsystem,  a  water 
management  subsystem,  a  thermal  management  subsystem  and 
a  power  electronics  subsystem  [1].  In  the  air  supply  subsystem, 
a  blower  or  a  compressor  is  generally  used  to  feed  air  to  the  cathode 
channel  of  a  stack.  As  the  supplied  air  travels  along  the  channel, 
oxygen  is  simultaneously  being  consumed  by  the  electrochemical 
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reaction.  The  oxygen  concentration  at  the  wall  catalyst  layer  is 
lower  than  that  in  the  core  cathode  flow  channel  due  to  the  electric 
current  harvested  at  the  wall  reaction  site.  The  fuel  cell  perfor¬ 
mance  is  strongly  affected  by  the  diffusive  mass  transport  through 
the  electrode  layer  which  is  proportional  to  the  reactant  concen¬ 
tration  gradient  between  the  catalyst  layer  in  the  membrane- 
electrode  assembly  (MEA)  and  the  flow  channel  in  the  bipolar  plate. 

As  the  current  generated  from  the  cell  increases,  the  higher 
electrochemical  reaction  rate  is  required.  Therefore,  the  reactant 
concentration  at  the  reaction  site  decreases  since  the  reactant  is 
rapidly  consumed  to  meet  the  increased  current  output.  The  mass 
transport  limit  is  reached  when  the  reactant  concentration  at  the 
catalyst  layer  drops  to  zero  [1,2].  Such  a  mass  transport  limit 
determines  the  limiting  current  density  ji  which  is  the  maximum 
current  output  from  the  cell.  Fuel  cells  cannot  operate  beyond  ji,  at 
which  the  reversible  cell  voltage  is  entirely  canceled  out  by  the 
concentration  loss.  In  the  fuel  cell  mass  transport  design,  therefore, 
the  focus  is  on  increasing  ji  which  is  determined  by  the  maximum 
reactant  transport  rate. 

To  extend  jL,  it  is  important  to  maintain  the  reactant  concen¬ 
tration  in  the  flow  channel  as  high  as  possible.  Thus,  high-powered 
blowers  for  enhanced  reactant  feeding  are  adopted  in  practical  PEM 
fuel  cell  systems.  However,  the  increase  of  blowing  capacity  brings 
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Nomenclature 

Sc 

Schmidt  number 

t 

time  (s) 

a 

pulsating  amplitude  (m) 

u 

velocity  in  x  direction  (m  s-1) 

A 

dimensionless  pulsating  amplitude 

U 

dimensionless  velocity  in  X  direction 

Ar 

aspect  ratio 

V 

voltage  (V) 

b 

half  of  cathode  channel  width  (m) 

W 

cathode  channel  width  (m) 

c 

oxygen  concentration  (mol  m-3) 

Wo 

Womersley  number 

C 

dimensionless  oxygen  concentration 

*>y 

axial  and  transverse  coordinates  (m) 

De 

Dm 

effective  time-averaged  dispersion  coefficient  (m2  s-1) 
oxygen  molecular  diffusivity  (m2  s-1) 

X,  Y 

dimensionless  axial  and  transverse  coordinates 

Ed 

dispersion  coefficient  ratio 

Subscripts  and  superscripts 

Ej 

limiting  current  density  enhancement  factor 

avg 

average 

EP 

maximum  power  density  enhancement  factor 

b 

bulk  state 

Ea 

wall  shear  stress  ratio 

cell 

fuel  cell 

f 

pulsating  frequency  (Hz) 

in 

channel  inlet 

F 

Faraday  constant  (96,485.34  C  mol-1) 

max 

maximum 

H 

cathode  channel  height  (m) 

out 

channel  outlet 

j 

current  density  (A  cm-2) 

St 

standard  condition,  273  K  and  1  atm 

k 

limiting  current  density  (A  cm-2) 

stack 

fuel  cell  stack 

k0 

pressure  gradient  in  steady  flow  (kg  m-2  s-2) 

0 

steady  component 

k\ 

l 

amplitude  of  pressure  fluctuation  (kgm-2  s-2) 
channel  length  occupied  by  MEA  (m) 

1 

unsteady  component 

L 

cathode  channel  length  (m) 

Greek 

symbols 

n 

number  of  mole  electrons  per  mole  reactant 

e 

amplitude  of  velocity  fluctuation 

N" 

axial  mass  flux  (mol  m-2  s-1) 

V 

efficiency 

P 

total  pressure  (kgm-1  s-2) 

P 

fuel  utilization  coefficient 

P 

power  density  (W  cm-2) 

V 

kinematic  viscosity  (m2  s-1) 

Qw 

oxygen  mass  flux  at  a  permeable  side  wall 

P 

fluid  density  (kgm-3) 

(mol  m-2  s-1) 

o 

wall  shear  stress  (kgm-1  s-2) 

Re 

Reynolds  number 

T 

dimensionless  time 

S 

interaction  function  between  velocity  and 
concentration  profiles 

(O 

angular  velocity  of  pulsating  flow  (s-1) 

forth  the  excessive  parasitic  power  loss  and  increased  system 
volume,  weight  and  noise  level.  For  this  reason,  a  new  mass 
transport  enhancing  scheme  has  been  suggested  to  reduce  high 
reactant  stoichiometry  [5],  and  they  demonstrated  that  the  fuel  cell 
performance  could  be  enhanced  by  the  reactant  pulsating  supply 
which  enhanced  axial  dispersion  of  the  reactant. 

This  paper  intends  to  address  the  impact  of  the  cathode  flow 
pulsation  on  the  PEM  fuel  cell  performance.  The  polarization  and 
corresponding  power  curves  are  experimentally  obtained  to  iden¬ 
tify  the  influence  of  the  Womersley  number  and  the  Reynolds 
number  in  the  cathode  flow  channel  on  the  overall  performance  of 
a  10-cell  PEM  fuel  cell  stack.  The  fuel  cell  performance  in  the  case 
of  a  pulsating  cathode  flow  is  compared  with  that  in  the  steady 
non-pulsating  flow  in  terms  of  the  maximum  power  density  and 
the  limiting  current  density.  In  addition,  a  theoretical  mass  trans¬ 
port  analysis  for  the  pulsating  cathode  flow  is  conducted  to  verify 
the  present  experimental  results  for  PEM  fuel  cell  stack. 

2.  Experimental  setup  and  procedure 

The  experimental  setup  was  composed  of  an  air  supply, 
a  hydrogen  supply,  a  nitrogen  supply,  a  temperature  control  unit 
and  a  PEM  fuel  cell  stack,  as  shown  in  Fig.  1.  An  NP50  PEM  fuel  cell 
stack  (heliocentric  Energiesysteme  GmbH)  with  10  cells  was 
employed.  The  active  area  of  a  cell  was  5  cm  x  5  cm  and  each  cell 
had  14  parallel  and  straight  open-air  cathode  channels  of  1.93  mm 
(H)  x  3.0  mm  ( W)  x  70  mm  (L)  as  sketched  in  Fig.  2. 

Compressed  hydrogen  of  99.999%  purity  was  supplied  to  the 
anode,  which  was  controlled  by  an  MFC.  Air  was  fed  to  the  cathode 


channels  from  a  compressed  air  tank  through  a  rectangular  Plexi¬ 
glas  duct  of  60  mm  in  width,  56  mm  in  height  and  870  mm  in  length, 
where  a  pulsating  component  was  added  by  a  woofer.  Three-stage 
mesh  screens  and  a  honeycomb  were  installed  to  produce 
a  uniform  inlet  flow  with  low  turbulence  intensity  of  0.68%  for  the 
steady  flow  at  Re  =  135.  The  flow  rate,  temperature  and  humidity 
were  controlled  by  an  MFC,  a  heater  and  a  humidifier,  respectively. 

An  acoustic  woofer  (Sammi  Sound  Technology  Corporation, 
SR-08B100)  was  installed  in  a  rectangular  chamber  to  produce  an 
oscillatory  flow.  The  rectangular  chamber  for  the  woofer  was  con¬ 
nected  to  the  Plexiglas  duct  with  a  flexible  tube  to  eliminate  trans¬ 
mission  of  structural  vibration.  A  function  generator  (HP-34401A) 
produced  a  sinusoidal  signal  with  a  specific  frequency  and  the  signal 
was  amplified  by  a  signal  amplifier  (Inkel,  AX7030G).  The  function 
generator  and  the  signal  amplifier  consumed  the  electric  power  up 
to  23  W  and  35  W,  respectively.  The  amplified  signal  was  delivered 
to  the  woofer  via  a  digital  oscilloscope  (LeCroy,  LT342)  to  check  the 
input  frequency  and  voltage.  The  time-variant  pulsating  velocity 
was  measured  at  the  mid  portion  of  the  cathode  channel  inlet  by 
a  hot-wire  anemometry  (Dantec  Streamline®  System)  to  determine 
the  amplitude  and  frequency  of  the  pulsating  flow  and  the  time- 
averaged  flow  velocity.  After  reaching  a  time-periodic  quasi 
steady-state,  the  experimental  data  were  stored  in  a  data  acquisition 
device  (Yokogawa,  DS  400)  for  analysis. 

In  an  effort  to  get  consistent  experimental  data,  the  following 
procedure  was  adopted  during  the  start-up,  performance 
measurement  and  shutdown  of  stack  operation  [5,6].  At  the 
beginning  of  the  experiment,  nitrogen  was  fed  to  the  stack  for 
20  min  to  purge  impurities  inside  the  anode  and  cathode  channels. 
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Fig.  1.  Schematic  diagram  of  the  experimental  setup. 


Then,  air  and  hydrogen  were  supplied  for  30  min  under  a  current 
loading  of  3  A  to  activate  the  MEA  inside  the  stack.  After  reaching 
a  steady-state,  the  open-circuit-voltage  was  measured.  The  current 
loading  was  increased  to  a  specific  value  by  using  an  electronic 
loader  (Daegil  Electronics  Co.,  EL-500P).  The  voltage  output  was 
determined  again  from  the  stack  in  the  steady-state  based  on  the 
galvanostatic  techniques  [1].  The  polarization  curves  were  finally 
obtained  for  various  current  loading  conditions.  After  finishing  the 
measurement,  the  external  loading  was  removed  and  the  air  and 
hydrogen  supplies  were  shut  off.  Then  the  stack  was  purged  with 
nitrogen  again  to  prevent  degradation  of  MEA  performance  due  to 
crossover  of  hydrogen  and  oxygen  remaining  at  the  anode  and 
cathode  channels. 


3.  Theoretical  analysis 


3.1.  Mathematical  formulation 


Consider  a  cathode  channel  (width  2b  and  length  l)  with  a  mass 
permeable  side  wall,  as  sketched  in  Fig.  3.  The  governing  equations 
for  a  one-dimensional,  unsteady,  fully  developed,  laminar  flow  of 
an  incompressible  fluid,  are  as  follows  [7,8] 


du  _  1  dp  6 2u 
dt  ~  ~p  9x  +  ?,0y2  ’ 


Oc  Oc  _ 

B  +  “  S"D 


9  2C 

ma72 


(1) 

(2) 


Fig.  2.  Configuration  of  the  cathode  channels. 


where  u{y,t)  represents  the  fluid  velocity  in  the  x  direction  and  c(x,y,t) 
denotes  the  oxygen  concentration  in  the  cathode  channel.  The  flow 
pulsation  is  driven  by  a  sinusoidal  pressure  pulse,  namely  [8-15] 

<3» 

where  k o  is  the  pressure  gradient  in  steady  flow,  k\  is  the  amplitude 
of  the  pressure  fluctuation,  i  denotes  and  the  subscript  R  refers 

to  the  real  part  of  the  terms  shown.  To  evaluate  k\,  it  is  assumed  that 
the  propagation  of  the  pulse  wave  in  the  channel  follows  the  same 
mechanism  as  the  propagation  of  a  sinusoidal  pressure  wave  in 
a  fluid,  namely  [9,10] 

/<i  =  (o2pa  ,  (4) 

where  w  is  the  angular  velocity  of  the  pulsating  flow,  a  is  the 
pulsating  amplitude,  and  p  is  the  fluid  density.  The  associated 
boundary  conditions  are  expressed  as  [1,2,8,11-13] 

=  Qwi  c0blx=0  =  cin- 

y=-b 

(5) 

Then,  Eqs.  (1)— (5)  are  non-dimensionalized  in  the  following 
fashion  [10-12] 


du 

ay 


=  °, 

=h  9y 
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Fig.  3.  Control  volume  for  a  pulsating  cathode  channel  flow. 


u  = 


uOavg 


X  = 


Dm* 

w0avg^2 

Re  .  M 


y  =  ;  -  -.Ji,  C  -  c  c‘n  , 

b  qwb/Dm 


OavgL 


cc  _  ^  A  —  — 
Dm’ 


(6) 


where  Wo,  Re,  Sc  and  A  are  the  Womersley  number,  the  Reynolds 
number,  the  Schmidt  number  and  the  dimensionless  pul¬ 
sating  amplitude,  respectively.  The  Womersley  number  is  the 
non-dimensional  oscillating  frequency  of  the  cathode  inlet  flow.  The 
Reynolds  number  denotes  the  non-dimensional  and  time-averaged 
air  flow  velocity  and  the  Schmidt  number  represents  the  relative 
magnitudes  of  molecular  momentum  diffusion  and  mass  diffusion 

rb 

in  the  flow.  Here,  u0aVg  =  (1/b)  /  uo  dy  =  (k0b2)/(3pr)  is  the 

cross-sectional  average  velocity  in  a  steady  flow,  and  Dm  is  the  bulk 
oxygen  molecular  diffusivity. 

The  governing  equations  and  the  associated  boundary  condi¬ 
tions,  in  properly  non-dimensionalized  form,  are 

dU-  3  2Wo^rir1  J_*v 

6t  Wo 2  Re  L  Jr  Wo2  8Y2  ’  v  ’ 


8C  1  a C  _  1  92C 

9t  +  Wo^Sc  OX  Wo2Sc  9Y2  ’ 


(8) 


—  1>  Q)blx=0  —  0, 

y=-i 

(9) 

where  Cob  is  the  bulk  concentration  in  steady  flow  defined  as  Q)b  - 
j\c0dY/  j\dY  [11-13], 

3.2.  Mass  transport  analysis 

Assuming  that  both  the  dimensionless  velocity  U  and  the 
dimensionless  concentration  C  consist  of  a  steady  component  and  an 
imposed  unsteady  component  [11-13],  i.e., 

U(Y,  t)  =  U0(Y)  +  U,(Y,z),  (10) 


a u 
ay 


=  0,  U\Y=±1=0,  g 
y=o  aY 


=  o  — 
5  ay 


c(x,y,T)  =  c0(x,y)  +  c1(y,i),  (li) 

and  substituting  Eqs.  (10)  and  (11)  into  Eqs.  (7)— (9),  the  analytic 
solutions  for  the  velocity  and  the  concentration  can  be  obtained 

U0(Y)  =  ~(y2-  l),  (12) 


2  Wo2A 
Re 


(13) 


C0(X,Y) 


X  Y^  _  3Y2  Y  39 
‘2+T6__8_+2  +  560’ 


(14) 


«'•')  -  ®KL' 

where 

R(Y)  =  1  fi  _  cosh(i1/2WoY)~ 

^  '  i  cosh  (i1/2  Wo) 

tanh  (i1 /2  wo)  cosh  (i1 /2  yyoSc1 /2  y) 
(1 -Sc)5c1/2sinh(i1/2WoSc1/2) 

cosh  (i1 /2  Wo  y)  1 

(1  -  5c)  cosh  (i1/2  Wo)  5c ' 


(15) 


(16) 


(17) 


Now,  it  is  possible  to  calculate  the  enhanced  mass  transport  in  the 
pulsating  flow  channel  from  the  velocity  and  concentration  profiles 
obtained.  The  effective  time-averaged  dispersion  coefficient  De  can 
be  determined  as  follows  [14,15] 


b 

-De^  +  T  J  u0(y)c0(x,y)  dy 

-b 

2tu/w  b 

=  -DmS  +  ^  /  ^/«(y.twx,y,t)dydt.  (18) 

0  -b 

The  left-hand  side  of  Eq.  (18)  denotes  the  effective  diffusion  flux  by 
the  flow  oscillation  plus  the  mass  flux  by  the  mean  flow,  and  the 
right-hand  side  represents  the  molecular  diffusion  flux  plus  the 
time-averaged  convective  mass  flux  produced  by  the  interaction  of 
the  cross-stream-varying  velocity  and  concentration  profiles. 
Substituting  Eqs.  (10)— (15)  into  Eq.  (18)  and  conducting  the  time- 
integration  yields  the  dispersion  coefficient  ratio  which  indicates 
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a  measure  of  the  effective  time-averaged  dispersion  coefficient  by 
flow  pulsation  in  comparison  to  the  bulk  molecular  diffusivity: 

c  De  1  ,  Sc2  Wo2/!2  f  (or,  d  V  W  }  Sc2Wo2A2 

£D=Dm=1+ — s —  J  yBG  +B  G)AY=J  =1+ — 8 — s> 

-1  -1 

(19) 

where  the  superscript  *  is  the  complex  conjugate  of  the  function. 
Here,  the  parameter  5  represents  the  interaction  of  the  velocity  and 
concentration  profiles. 

To  assess  how  much  the  flow  pulsation  augments  the  mass 
transport  in  the  channel,  an  analogy  between  mass  and  momentum 
transfer  is  implemented.  It  has  been  well  reported  that  the  velocity 
profiles  exhibit  the  Hagen— Poiseuille  shape  at  a  low  Wo.  However, 
the  parabolic  shape  is  lost  and  the  thinner  boundary  layer  is  shown, 
as  Wo  increases  [8,10].  Thus,  the  augmented  velocity  gradient  in  the 
near-wall  region  brings  forth  a  beneficial  effect  on  mass  transfer 
[9,10].  The  time-averaged  shear  stress  at  the  wall  in  the  pulsating 
flow  can  be  written  as  [9] 

2tz/(x) 

°  =  /  (S)(uavgW)dt,  (20) 

o 

where  uavg(t)  denotes  the  cross-sectional  average  velocity  at  any 
time,  and  uavg  is  the  time-averaged  bulk  mean  velocity,  namely 

b  2tu/w 

Wavg(t)  —  g  J  u(y,t)dy  ,  Uavg  =  ^  J  tfaVg(t)d t.  (21) 

0  0 


Substituting  Eqs.  (3)  and  (10)  into  Eq.  (20),  the  time-averaged 
wall  shear  stress  ratio  of  the  pulsating  flow  to  the  steady  flow  can 
be  drawn  as  follows: 


a  2 Wo5 A2  sinh(V2Wo)  -  sin(V2  Wo) 

0b  +  3V2 Re2  cosh  (y/2' Wo)  +  cos  (VZ Wo) 


(22) 


Through  an  analogy  between  mass  and  momentum  transfer,  the 
limiting  current  density  enhancement  factor  which  means  the 
mass  transport  ratio  of  the  pulsating  flow  to  the  steady  flow  can  be 
expressed  in  terms  of  the  time-averaged  wall  shear  stress  as  follows 
[9,10]: 


n_=  (r\V3=  h. 

Nq  VW  JoL 


(23) 


It  is  clearly  discernible  from  Eqs.  (22)  and  (23)  that  the  limiting 
current  enhancement  by  flow  pulsation  is  pronounced  as  the 
Womersley  number  Wo  and  the  pulsating  amplitude  A  increase. 
The  flow  pulsation  effect,  on  the  other  hand,  dwindles  as  the  Rey¬ 
nolds  number  Re  increases. 

In  the  present  mass  transport  analysis,  the  parameter  ranges  are 
50  <Re<  150,  0.1  <  Wo  <  10, 1  <A  <  1  and  Sc  =  0.765,  which  cover 
the  parameter  ranges  in  the  present  fuel  cell  experiments.  It  should 
be  noted  that  the  tidal  displacement  should  be  smaller  than  the 
flow  channel  length  to  prevent  direct  convective  mass  transport  in 
the  axial  direction  [15,16],  and  this  limitation  can  be  expressed  as 
A  <  16.67  for  the  present  pulsating  cathode  channel. 


4.  Results  and  discussion 


4.1.  Velocity  measurement 


represents  the  time-averaged  inlet  velocity  which  is  related  to  the 
steady  mean  flow,  and  e  is  the  amplitude  of  the  velocity  fluctua¬ 
tion  related  to  the  flow  oscillation.  The  tidal  displacement  repre¬ 
sents  the  cross-stream-averaged  maximum  axial  distance  which 
the  fluid  elements  travel  during  a  one-half  period  of  the  pulsation 
[15].  The  pulsating  amplitude  is  half  of  the  tidal  displacement  [16], 
namely 


To  characterize  the  pulsating  flow,  it  is  useful  to  introduce  the 
time-variant  pulsating  velocity  u  =  u0  +  esin(wt).  Here,  Uo 


Fig.  4.  Time-variant  pulsating  velocity  u  (a  =  0.0025  m):  (a )/=  10  Hz,  (b)/=  30  Hz  and 
(c)/=  50  Hz. 
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a  =  -  J  £*sin(wt)dt  =  e/o). 
o 


(24) 


The  Womersley  number  Wo  =  h(w/^)^2  and  the  well-known 
Reynolds  number  Re  =  2u0avgb/v  can  be  calculated  from  the 
measured  velocity.  Here,  a  polynomial  approximation  for  the 
maximum-to-average  velocity  ratio  in  laminar  flow  through 
a  rectangular  duct  is  adopted  to  predict  the  cross-stream  averaged 
velocity  Uoavg  by  using  the  measured  center-line  velocity  [17], 


=  1  (l  +  0.546688Ar  + 1 .55201 3Ar2  -  4.059427Ar3 

w0avg  ^  ' 

+  3.214927Ar4-0.857313Ar5) .  (25) 


The  measured  time-variant  pulsating  velocities  at  the  mid  portion 
of  the  cathode  channel  inlet  are  displayed  in  Fig.  4  when  the 
pulsating  amplitude  is  fixed  at  a  =  0.0025  m.  The  time-averaged 
inlet  velocity  uo  is  varied  from  0.87  m  s-1  to  1.53  m  s-1  and  the 
amplitude  of  pulsating  velocity  e  ranges  from  0.16 ms-1  to 
0.78  m  s-1  to  fix  the  pulsating  amplitude  for  all  the  flow  conditions 
tested.  The  frequency  spectrum  for  the  velocity  fluctuation  at 
Re  =  112  discloses  that  the  present  acoustic  excitation  can  provide 
quite  accurate  sinusoidal  flow  pulsations  at  a  specific  forcing 
frequency /and  pulsating  amplitude  a,  see  Fig.  5. 


4.2.  Fuel  cell  performance 

The  polarization  and  corresponding  power  curves  for  the  NP50 
PEM  stack  according  to  the  Womersley  number  Wo  are  demon¬ 
strated  in  Fig.  6.  In  order  to  scrutinize  the  effect  of  Wo,  the  pulsating 
amplitude  is  maintained  at  a  =  0.0025  m  and  the  inlet  gas 
temperatures  at  the  anode  and  the  cathode  are  fixed  at  30  °C  and 
35  °C,  respectively.  Also,  the  relative  humidities  at  the  anode  and 
cathode  inlet  are  set  at  70%  and  80%,  respectively.  Such  an  experi¬ 
mental  condition  may  be  unrealistic  comparing  to  the  actual  air  and 
fuel  supply  devices  in  PEM  fuel  cell  systems.  However,  the  impact  of 
the  flow  parameters  as  well  as  the  oscillating  frequency  on  the 
performance  of  the  PEM  fuel  cell  can  be  clearly  and  quantitatively 
identified  by  adopting  the  present  experimental  setup.  The  Rey¬ 
nolds  number  of  cathode  flow  is  changed  in  the  range  of  77  (Fig.  6a) 
to  135  (Fig.  6c)  at  a  constant  hydrogen  flow  rate  of  2  1pm  and  a  stack 
temperature  of  50  °C. 


Fig.  5.  Power  spectrum  of  the  forcing  frequency  (Re  =  112). 
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Fig.  6.  Polarization  and  power  curves  according  to  the  Womersley  number:  (a)  Re  =  77, 
(b)  Re=  112  and  (c)  Re=  135. 


It  is  evident  that  the  enhancement  of  Vstack  by  the  flow  pulsation 
is  negligible  at  low  current  loadings.  As  the  Womersley  number 
increases,  however,  the  cathode  flow  pulsation  significantly 
enhances  the  fuel  cell  performance  at  the  concentration  loss  region. 
The  concentration  loss  in  the  polarization  behavior  is  dramatically 
reduced  by  the  cathode  flow  pulsation,  while  the  activation  and 
ohmic  losses  remain  unchanged.  The  limiting  current  density  ji 
increases  with  Wo  due  to  the  enhanced  dispersion  of  reactant  in  the 
cathode  channels.  Accordingly,  the  maximum  power  density  is  also 
augmented.  As  shown  in  Fig.  6a,  the  limiting  current  density  and 
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the  maximum  power  density  for  the  steady  flow  at  Re  =  77  is 
0.4  A  cm-2  and  1.62  W  cm-2,  respectively.  As  the  oscillating  flow  is 
added  at  Wo  =  6.55  and  a  =  0.0025  m,  however,  they  are  dramati¬ 
cally  augmented  up  to  0.56  A  cm-2  and  2.20  W  cm-2,  respectively. 
The  effect  of  cathode  flow  pulsation  at  Re  =  112  is  given  in  Fig.  6b. 
The  limiting  current  and  maximum  power  densities  are 
0.60  A  cm-2  and  2.26  Won-2  for  a  steady  non-pulsating  flow, 
respectively.  They  are  increased  up  to  0.74  A  cm-2  and  2.54  W  cm-2 
by  the  cathode  flow  pulsation  of  Wo  =  6.55,  respectively.  For 
Re  =  135,  the  limiting  current  densities  are  0.74  A  cm”2  and 
0.86  A  cm-2  for  steady  flow  and  pulsating  flow  of  Wo  =  6.55  oper¬ 
ations,  respectively,  see  Fig.  6c.  The  maximum  power  density  is 
augmented  from  2.53  W  cm-2  to  2.72  W  cm”2  It  is  clearly 
discernible  that  the  enhancement  of  the  limiting  current  density 
and  the  maximum  power  density  dwindles  away  as  Re  increases. 

To  assess  the  effect  of  the  cathode  flow  pulsation  on  the  overall 
fuel  cell  performance,  it  is  convenient  to  introduce  the  enhance¬ 
ment  factors  in  the  limiting  current  density  and  the  maximum 
power  density  as  follows: 


=  Jl/Jol  j 

(26) 

Ep  =  Emax/Pomax  ? 

(27) 

where  the  subscript  0  denotes  the  steady  non-pulsating  flow.  Thus, 
the  enhancement  factor  represents  the  ratio  of  the  pulsating 
cathode  flow  to  the  steady  flow.  Fig.  7  illustrates  the  variation  of  the 
limiting  current  density  enhancement  factor  Ej  and  the  maximum 
power  density  enhancement  factor  Ep.  The  enhancement  factors 
increase  as  the  cathode  inlet  flow  oscillates  at  higher  Wo.  The 
pulsating  cathode  flow  of  Wo  =  6.55  gives  rise  to  40.0%,  23.3%  and 
16.2%  increments  in  the  limiting  current  density  when  the  Reynolds 
numbers  are  77,  112  and  135,  respectively,  see  Fig.  7a.  As  high¬ 
lighted  in  Fig.  7b,  it  also  provides  35.5%,  12.5%  and  7.9%  higher 
maximum  power  density  compared  to  the  steady  non-pulsating 
flow  for  Re  =  77,  112  and  135,  respectively.  As  expected,  the 
enhancement  factors  are  more  pronounced  as  Re  decreases.  These 
results  imply  that  the  flow  pulsation  is  more  effective  at  low  Re.  For 
high  Re,  the  oxygen  concentration  in  the  cathode  channel  is  suffi¬ 
cient  for  the  full  electrochemical  reaction  so  that  the  enhancement 
by  the  flow  pulsation  is  meager.  It  is  anticipated  that  the  trend  of 
enhanced  performance  by  the  cathode  flow  pulsation  can  be  also 
obtained  in  actual  PEM  fuel  cell  systems  by  adopting  an  appropriate 
pulse  generating  device. 

4.3.  Fuel  cell  efficiency 

The  fuel  cell  efficiency  based  on  the  higher  heating  value  is 
expressed  by  [3]: 


where,  Vcen  is  the  average  cell  voltage  and  p  is  the  fuel  utilization 
coefficient  defined  as  the  fuel  ratio  of  the  consumption  in  a  cell  to 
the  supply  into  a  cell.  The  efficiency  of  a  fuel  cell  can  be  estimated 
from  the  simple  measurement  of  VCeii  when  a  good  estimate  for  g  is 
0.95  [3]. 

The  fuel  cell  efficiency  p  as  a  function  of  power  density  P  at 
a  pulsating  amplitude  of  a  =  0.0025  m  is  presented  in  Fig.  8.  Fig.  8a 
shows  that  a  maximum  power  density  of  1.62  W  cm-2  is  obtained  at 
the  efficiency  of  37.2%  for  the  steady  non-pulsating  flow  of  Re  =  77. 
The  maximum  power  density  increases  up  to  1.74  W  cm-2, 
1.98  W  cm-2  and  2.20  W  cm-2  for  the  Womersley  number  of  2.93, 
5.07  and  6.55,  respectively.  The  corresponding  fuel  cell  efficiencies 
are  34.8%,  31.7%  and  29.4%,  respectively.  Therefore,  the  substantially 


Fig.  7.  Enhancement  factors  versus  the  Womersley  number:  (a)  limiting  current 
density  and  (b)  maximum  power  density. 

enhanced  maximum  power  density  with  Wo  can  be  achieved  at 
a  slightly  reduced  efficiency  by  adopting  the  cathode  flow  pulsation. 
Even  at  the  identical  efficiency,  the  cathode  flow  pulsation  yields  the 
prominent  enhancement  of  power  density.  It  should  be  noticed  in 
Fig.  8b  and  c  that  the  maximum  power  density  enhancement  with 
Wo  dwindles  as  the  Reynolds  number  increases.  Flence,  it  is  clear 
that  the  cathode  flow  pulsation  is  more  effective  at  low  Re. 

4.4.  Mass  transport  in  a  pulsating  cathode  flow  channel 

The  dispersion  coefficient  ratio  can  be  estimated  from  the 
experimental  results  of  the  PEM  fuel  cell  stack  with  a  cathode 
pulsating  flow.  Assuming  the  time-averaged  axial  concentration 
gradient  is  constant,  the  time-averaged  axial  mass  flux  can  be 
expressed  as 

K  =  U0avg(C0in  -  C0out)  -  DmC°out  ~  C°in  =  £ ,  (29) 

U  =  Uavg(Cjn  —  Cout)  —  De- - y  =  -j-p,  (30) 

for  the  steady  non-pulsating  flow  and  the  pulsating  flow,  respec¬ 
tively.  Flere,  l  is  the  flow  channel  length  occupied  by  the  MEA,  n  is 
number  of  electrons  transferred  in  the  reaction,  F  is  Faraday’s 
constant  (F=  96485.34  C  mol-1),  and  the  bars  designate  the  time- 
averaged  value.  Supposing  that  all  the  oxygen  in  the  cathode 
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Fig.  8.  Fuel  cell  efficiency  according  to  the  Womersley  number:  (a)  Re  =  77, 
(b)  Re  =  112  and  (c)  Re=  135. 


channel  is  entirely  consumed  by  the  electrochemical  reaction  in  the 
mass  transport  limit  situation,  the  time-averaged  bulk  concentra¬ 
tion  at  the  channel  outlet  can  be  assumed  to  be  zero  [1,2].  Sub¬ 
tracting  Eq.  (30)  from  Eq.  (29)  gives  the  dispersion  coefficient  ratio 
with  the  equalities  of  u0aVg  =  Uavg  and  c0in  =  cin: 


De  =  i  +  l  Jl-Jol 
Dm  Coin  nFDm 


(31) 


The  geometrical  parameter  and  the  physical  properties  are  listed  in 
Table  1.  Here,  the  oxygen  molecular  diffusivity  Dm  is  a  function  of 
temperature  and  pressure  [2,7]. 


Table  1 

Geometrical  parameters  and  physical  properties. 


Quantity 

Value 

Reference 

Cathode  channel  length 

0.05  m 

Present 

occupied  by  MEA,  l 

experiment 

Cathode  channel  width,  W  =  2b 

0.003  m 

Pulsating  amplitude,  a 

0.0025  m 

Dimensionless  pulsating 

1.65 

amplitude,  A  =  a/b 

Cathode  pressure,  p 

1  atm 

Cathode  channel  temperature,  T 

35  °C 

Relative  humidity  of  inlet  air,  RH 

80% 

Number  of  mole  electrons  transferred 

4 

[1-3] 

in  the  cathode  reaction,  n 

Oxygen  diffusivity  in  gas  at 

0.18  x  10"4  m2  s— 1 

[2] 

standard  condition,  Dst 

Oxygen  molecular  diffusivity, 

2.157  x  10"5  m2  s— 1 

[2,7] 

Dm  =  Dst(T/Tst)3/2(Pst/p) 

Oxygen  concentration  at 

7.945  mol  m  3 

[18] 

the  channel  inlet,  c0in 

Kinematic  viscosity  of  inlet  air,  v 

1.650  x  10  5  m2  s  1 

[18] 

Schmidt  number,  Sc  =  v/Dm 

0.765 

[2,7,18] 

Fig.  10.  Relation  between  the  limiting  current  density  and  the  time-averaged  wall 
shear  stress. 
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In  Fig.  9,  the  dispersion  coefficient  ratios  evaluated  from  the 
fuel  cell  experiments  [Eq.  (31)]  are  compared  with  the  theoretical 
analysis  [Eq.  (19)]  when  the  dimensionless  pulsating  amplitude 
and  the  Schmidt  number  are  A  =  1.65  and  Sc  =  0.765,  respectively. 
The  dispersion  coefficient  ratio  increases  with  the  increase  of  Wo. 


Y 


Y 
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Fig.  11.  Effect  of  the  Womersley  number  on  the  time- variant  velocity  profile  (A  =  7  and 
Re  =  50):  (a)  Wo  =  0.1,  (b)  Wo  =  1  and  (c)  Wo  =  10. 


The  experimental  results  show  that  the  effect  of  Wo  is  in  well 
accordance  with  the  theoretical  analysis.  Some  discrepancy  is  due 
to  the  fact  that  the  present  analysis  adopts  the  simple  constant 
mass  flux  boundary  condition,  excludes  the  reaction  site  where  the 
complicated  electrochemical  reaction  occurs,  and  disregards  the 
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Fig.  12.  Effect  of  the  Womersley  number  on  the  time-variant  concentration  profile  at 
the  channel  inlet  (A  =  7  and  Re  =  50):  (a)  Wo  =  0.1,  (b)  Wo  =  1  and  (c)  Wo  =  10. 
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3D  effect  of  the  diffusive  transport  into  the  catalyst  layer,  to  name 
a  few. 

In  Fig.  10,  the  experimental  values  of  the  limiting  current 
density  enhancement  factor  Ej  are  plotted  against  the  time- 
averaged  wall  shear  stress  ratio  to  the  power  of  one-third  E ]/3 
evaluated  from  Eq.  (22).  Also,  the  theoretical  prediction  given  by 
Eq.  (23)  is  presented  as  a  solid  line.  The  experimental  results  are 
shown  to  be  consistent  with  the  theoretical  ones  within  8.6% 
deviation.  It  is  evident  that  flow  pulsation  contributes  to  increase 
the  limiting  current  density. 

The  analytical  velocity  profiles  according  to  Wo  during  a  single 
cycle  of  the  flow  pulsation  are  displayed  in  Fig.  11.  As  the 
Womersley  number  Wo  increases  from  0.1  to  10  at  A  =  7  and 
Re  =  50,  the  flow  field  is  substantially  affected  by  the  flow  pulsation 
[8,13].  The  velocity  magnitude  is  augmented  and  the  boundary 
layer  becomes  thinner  with  the  increase  of  Wo.  The  velocity  profile 
at  the  core  becomes  flat. 

Fig.  12  illustrates  the  effect  of  flow  pulsation  on  the  concentra¬ 
tion  profiles  at  the  channel  inlet  (X  =  0)  when  the  dimensionless 
pulsating  amplitude  and  the  Reynolds  number  are  A  =  7  and 
Re  =  50,  respectively.  The  left  side  wall  shows  the  lowest  oxygen 
concentration  because  the  oxygen  gas  is  being  consumed  by  the 
electrochemical  reaction,  as  depicted  in  Fig.  3.  In  the  same  manner 
as  in  the  velocity  profile,  the  effect  of  flow  pulsation  on  the 
concentration  field  is  meager  when  Wo  is  low.  As  Wo  increases, 


Fig.  13.  Effective  time-averaged  dispersion  coefficient:  (a)  interaction  function 
between  the  velocity  and  concentration  profiles  S  and  (b)  dispersion  coefficient  ratio  ED. 


however,  the  concentration  profile  begins  to  oscillate  periodically 
because  the  pulsating  velocity  directly  affects  the  concentration 
field,  see  Eq.  (2).  It  is  discernible  from  Fig.  12c  that  the  concentra¬ 
tion  intensively  fluctuates  in  the  core  region  where  the  velocity 
oscillating  magnitude  is  high.  Also,  the  phase  lag  with  respect  to  the 
velocity  fluctuation  is  shown,  which  plays  a  key  role  in  mass 
transport  enhancement  by  the  interaction  of  the  velocity  and 
concentration  profiles  [14,15]. 


Fig.  14.  Limiting  current  density  enhancement  factor:  (a)  A  =  1,  (b)  A  =  4  and  (c)  A  =  7. 
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The  effective  time-averaged  dispersion  coefficient  is  elucidated 
from  Fig.  13  by  plotting  the  parameter  S  which  represents  the 
interaction  of  the  velocity  and  concentration  profiles.  The  interac¬ 
tion  function  5  increases  to  a  maximum  and  gradually  decreases 
beyond  Wo  ~  3  as  Wo  increases,  see  Fig.  13a.  As  shown  in  Fig.  13b, 
therefore,  the  increasing  rate  of  the  dispersion  coefficient  ratio  ED 
slows  in  the  region  of  Wo  =  1  because  the  dispersion  coefficient 
enhancement  is  proportional  to  the  product  of  the  Wo2  and  S,  see 
Eq.  (19).  Physically,  the  fluid  never  has  enough  time  for  diffusion 
from  the  high  concentration  core  region  to  the  low  concentration 
near-wall  region  or  vice  versa  under  fast  oscillation,  i.e.  high  Wo, 
which  means  the  molecular  diffusion  time  scale  is  much  larger  than 
the  one-half  period  of  the  oscillation  [15,19].  The  increase  of  the 
dimensionless  pulsating  amplitude  A  results  in  a  substantial 
enhancement  of  the  dispersion  coefficient  due  to  the  intensified 
pulsating  effect. 

Fig.  14  exhibits  the  limiting  current  density  enhancement 
factor  Ej  calculated  from  Eqs.  (22)  and  (23).  The  limiting  current 
density  enhancement  factor  is  almost  unity  which  is  equal  to  the 
steady  flow  value  at  low  Womersley  numbers.  As  Wo  increases, 
however,  it  is  steeply  increased  by  the  enhanced  axial  dispersion 
in  the  cathode  flow  channel.  The  enhancement  of  the  limiting 
current  density  is  more  pronounced  at  higher  A  because  of 
the  intensified  pulsating  effect.  Also,  it  is  discernible  that  the 
enhancement  factor  dwindles  as  Re  increases.  For  high  Re ,  the 
mass  transport  by  the  steady  forced  convection  in  the  channel  can 
be  enough  for  the  electrochemical  reaction  even  in  the  absence  of 
the  flow  pulsation. 

5.  Conclusions 

The  effect  of  the  pulsating  cathode  flow  on  the  overall  perfor¬ 
mance  of  a  PEM  fuel  cell  has  been  investigated.  The  polarization  and 
corresponding  power  curves  are  experimentally  obtained  to 
investigate  the  impact  of  the  Womersley  number  and  the  Reynolds 
number  on  the  performance  of  a  10-cell  PEM  fuel  cell  stack.  As  the 
pulsating  component  is  delivered  only  to  the  cathode  mainstream 
flow,  the  performance  of  the  PEM  fuel  cell  stack  in  the  concentration 
loss  region  is  prominently  increased  due  to  the  enhanced  axial 
dispersion.  The  increased  limiting  current  density  is  obtained  at 
higher  Wo.  The  cathode  flow  pulsation  at  Wo  =  6.55  yields  40.0%, 
23.3%  and  16.2%  augmented  limiting  current  densities  compared  to 
the  steady  non-pulsating  flow  for  Re  =  77, 112  and  135,  respectively. 
Accordingly,  the  maximum  power  density  increases  with  Wo  at 


a  slightly  reduced  fuel  cell  efficiency.  The  maximum  power  density 
shows  35.5%,  12.5%  and  7.9%  increments,  respectively.  Enhancement 
of  the  overall  performance  is  more  pronounced  at  low  Re.  Conse¬ 
quently,  the  flow  pulsation  extends  the  mass  transport  limit  in  the 
cathode  channel,  similarly  the  increase  of  steady  flow  rate  by 
adopting  high  capacity  blowers  or  compressors. 

In  addition,  a  theoretical  mass  transport  analysis  in  the  pulsating 
cathode  flow  channel  has  been  implemented  to  verify  the  experi¬ 
mental  results.  The  cathode  flow  pulsation  increases  the  oxygen 
axial  dispersion  in  the  cathode  channel  so  that  the  limiting  current 
density  increases  with  the  Womersley  number  and  the  pulsating 
amplitude.  As  the  Reynolds  number  increases,  however,  the 
enhancement  of  the  limiting  current  density  dwindles  due  to  the 
reduced  pulsating  effect.  This  theoretical  mass  transport  analysis 
shows  that  the  effect  of  the  relevant  parameters  is  in  well  accor¬ 
dance  with  the  10-cell  PEM  fuel  cell  stack  experiments. 
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